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ABSTRACT 

We derive self-similar solutions including cooling and heating in an Einstein de-Sitter 
^ . universe, and investigate the effects of cooling and heating on the gas density and tem- 

perature distributions. We assume that the cooling rate has a power-law dependence 
on the gas density and temperature, Aocp^T^, and the heating rate is TocpT. The 
values of A and B are chosen by requiring that the cooling time is proportional to the 
Hubble time in order to obtain similarity solutions. In the region where the cooling 
rate is greater than the heating rate, a cooling inflow is established, and the gas is 
^ ■ compressed and heats up. Because the compression is greater in the inner region than 

[ in the outer region, the temperature becomes an increasing profile toward the center. 

In particular, when a large infall velocity is produced due to an enormous energy loss, 
the slope of the density approaches a value that depends on A, B, and the velocity 
slope, and the slope of the temperature approaches —1. On the other hand, in the 
region where the heating rate is greater than the cooling rate, the infall velocity is 
suppressed, compression of the gas is weakened, and the gas cools down. The slope of 
the density becomes shallow due to suppression of the contraction, and the tempera- 
ture is lower than that without heating. The self-similar collapse presented here gives 
insights to the effects of cooling and heating on the gas distributions in galaxies and 
Q . clusters of galaxies. 

Key words: cosmology: theory - galaxies: clusters: general - gravitation - hydrody- 
namics - intergalactic medium - methods: analytical - radiation mechanisms: general. 
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"1 INTRODUCTION Sitter universe, Bertschinger(1985) derived the similarity so- 
^ lution of a top-hat density perturbation, n=3, where n is a 
Clusters of galaxies contain large quantities of hot gas. The power spectrum index, and showed that the density has a 
gas collapses under the influence of gravity, and loses ther- power- law profile against the radius, pocr"^-^^ Chuzhoy & 
mal energy due to radiative cooling. Especially, the effect Nusser(2000) obtained the asymptotic behaviour of similar- 
of cooling is important in the inner region where the cool- ity variables of adiabatic gas, and found that they strongly 
ing time is less than the Hubble time. Inside this region, depend on the initial power spectrum. For n >-2, they 
it is expected that the internal pressure decreases, the gas found that the density and temperature asymptotically ap- 
cools down and cooling flows occur. However, from recent proach pocr"'^'"'*'^''''"'^^-' and Tocr~^"~^^^^"'^^^ in the limit 
observations of clusters of galaxies, the expected cooled gas ^^0, respectively. For a self-similar solution with cooling, 
is not detected in the central region (Peterson et al. 2001; Abadi et al.(2000) obtained a similarity solution of coUi- 
Tamura et al. 2001). Furthermore, X-ray observations indi- gional gas in an Einstein de-Sitter universe, and used it to 
cate that the observed luminosity-temperature relation dif- evaluate the ability of SPH simulations. In their solution, 
fers from the self-similar relation (Arnaud & Evrard 1999; they assumed the cooUng function Axp^^^T. This power- 
Helsdon & Ponman 2000). As an origin of this discrepancy, law dependence is determined by the requirement that the 
it is supported that the gas is signiflcantly affected by non- cooling time must have a fixed fraction to the dynamical 
gravitational heating (Kaiser 1991; Evrard & Henry 1991). time in order to obtain the similarity solution. They found 
It is very important to investigate the influences of cooling that radiative cooling decreases the pressure support, and 
and heating on the gas distributions, and it would provide establishes a cooling inflow. As a result, the profiles of fiuid 
an understanding of the thermal evolution of the intracluster variables differ from those of the adiabatic solution. In par- 
medium, ticular, when the energy loss is large enough to produce a 
For adiabatic similarity solutions in an Einstein de- 
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large infall velocity, they found that the density and temper- 
ature approach p oc and T oc as r -^0, respectively. 

The purpose of this study is to investigate the effects of 
cooling and heating on the gas distributions using the simi- 
larity solution. The cooling rate per unit volume is assumed 
to be given by the power-law dependences on the density 
and temperature, Aocp^T^ , where A and B are related to 
the initial power spectrum. The heating rate is assumed to 
be proportional to the gas density and temperature, TocpT. 
Of course, the similarity solution is not directly applicable 
to the real structure, and cannot include detailed physical 
processes. However, it is possible to show exact structures 
of flows and to give clues concerning the effects of cooling 
and heating on the gas distributions of galaxies and clusters 
of galaxies. 

The layout of the paper is as follows. In Section 2, we 
present our model. Our results are shown in Section 3. Fi- 
nally, we discuss and summarize our results in Section 4. 



2 MODEL 

2.1 Spherical collapse in an expanding universe 

We consider here that an over-dense region has a density 
contrast (Si) at an initial Hubble time of tH,i, and that the 
background universe is an Einstein de-Sitter universe. The 
gravitational collapse in the flat universe becomes a simi- 
larity evolution because bounded shells extend to infinity. 
Fillmore & Goldreich(1984) introduced a scale-free initial 
density contrast, 



Mi ^mJ ^Ro' ' 



(1) 



where Mi — {4/3)'KpH,iri is the mass contained within the 
proper radius r^, pH,i = l/(67rGt|f j) is the background den- 
sity at tH,i, and Mo={4:/3)TvpH,iRo is a reference mass. Since 
5Mi must be constant or increase with increasing radius, and 
5i has to decrease with increasing radius, e is restricted to 
< e < 1. The variance of the density fluctuation is 

\5{x,tH)f ^ ^ aitH^M-"^, (2) 

where x is the comoving radius, and \5k\^ oc a^tuYk^ is the 
power spectrum, where a{tH) is an expansion factor. The 
relation between the index n and the parameter e is 



n = 3(2e - 1). 



(3) 



We assume that the initial velocity of matter is pure Hubble 
flow, and that the initial gas pressure is zero. The expand- 
ing matter slows down with the growth of the perturbation. 
Eventually, it reaches a maximum radius and decouples from 
the Hubble expansion. Since a(tir) oc , the perturbation 
grows as |5(r,tff)p oc a(te)"+V-("+3) oc tW)/3^-(n+3)^ 
based on equation (2). The density contrast at the turn- 
around radius is (5(na, iH)=(37r/4)^-l'-4.55(Peebles 1980). 
Thus, the time evolution of the turn-around radius and the 
enclosed mass are 



na(tif) oc t\j, 
2 

mta{tH) oc tj, 

where 



(4) 
(5) 



(6) 



After the turn-around, particles fall toward the center. The 
collapse of matter causes an increase of the central den- 
sity. Since the gas particles are decelerated, a shock wave 
is formed and propagates outward. The system undergoes 
self-similar evolution at a later time. 



2.2 Cooling and heating laws in similarity 

evolutions 

In the case of the collapse of adiabatic gas, the physical scale 
of the system is only determined by the gravitational dynam- 
ics. When radiative cooling and heating are included, the 
similarity evolution is broken due to the existence of many 
physical scales. In order to maintain the similarity evolution, 
it is required that all of the physical scales are proportional 
to the dynamical scale. The scale-free condition was con- 
structed by ensuring that the cooling time is proportional 
to the dynamical time (Owen et al. 1998; Abadi et al. 2000). 

We assume that the cooling rate per unit volume is 
A=Aop^C3^. Here, A and B are free parameters, Ao is the 
cooling coefficient, p is the gas density, and Cs is the sound 
speed. The cooling time becomes 

E A _ pi-^c^^^""®' 

= ^ =£;(-)- = ^^^-^, (7) 

where E is the specific thermal energy of the gas, and 7 = 
5/3 is the adiabatic index. The Hubble time evolves as 



o ao 



(8) 



where oo and Hq are the expansion factor and Hubble con- 
stant at the present epoch, respectively. The sound speed 

and density evolve as 



ao 

I ^ ^-3 
P=( — ) Po, 
ao 



(9) 
(10) 



where Cs,o and po are the sound speed and the gas density at 
the present epoch, respectively. We assume that the Hubble 
time is proportional to the cooling time by the constant tc, 
in order to maintain scale-free conditions (i.e., tH=tctcooi), 



Ih = t, 



p^-^cT-"^ 



(11) 



7(7 - l)Ao 

Substituting equations (8)- (10), the above equation becomes 



o „1-^^2{1-B) 

s""" ' 7(7-l)A<. ' 



(12) 



The a-dependence term must vanish in order to keep tc con- 
stant. Using equation (3), we find 



(1 - S) = 0. 



(13) 



Figure 1 displays the parameter space of A and B in 
relation to e. The shaded areas correspond to < e < 1. 
The solid lines display the lines of e=l, 2/3, 1/3, and 0. 
Abadi et al.(2000) selected at the point (^,B)=(3/2,1), 
which is independent of e. The points of free-free emission. 
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Figure 1. Relation between (A, B) and e. The shaded areas indi- 
cate the inside of the range < e < 1 . The solid lines are constant 
values of e=l, 2/3, 1/3, and 0. The closed circles represent the 
parameter positions that we calculated in this paper. The points 
of free-free emission and line emission are indicated by crosses. 



{A,B)={2,l/2), and line cooling, {A,B) = {2,~l/2), are indi- 
cated by crosses. 

Next, we assume that the heating rate per unit volume 
is given by 



■n 2 
r = ape, , 



(14) 



where a={tH,o/tH)ao, tH,o is the current Hubble time, and 
ao is the model parameter. The heating time becomes 



DE 
Dt 



(15) 



We assume that the heating time is proportional to the Hub- 
ble time by the constant th(i.e., tH=ththeat), 

I"" = 7(7 -Vo- 



2.3 Basic equations of collisional gas 

Neglecting the effect of dark matter, the basic equations 
with additional terms of cooling and heating are as follows: 

dm 2 
— = ATvr p, 
or 



dp 
dt 



dv dv 
dt dr 



d , 2n Gm 
^(pcj 5-, 



(17) 
(18) 
(19) 



1 del (?s dp . A-l 2B , 2 

— -77 =-7^0/5 Cs +10lCs, 



(20) 



7 — 1 dt p dt 

where v is the velocity, G is the gravitational constant, and 
m is the mass of the gas. 

To derive self-similar solutions, we introduce the simi- 
larity variable X and similarity functions, which are given 
by the following forms: 
r 

X = —, 

rta 

p = PhD, 



V = —V, 

tH 

Cs = — C, 

tn 

4 3 

m = -TVpHVtaM, 

n 

where pH=^/ {Qi^Gt^j) is the mean density of the universe, 
and Ti. is the same as aotnfi or 3i/oihiw,o/27(7 — 1) from 
equation (16). The density and the sound speed are ex- 
pressed as 

PH 



(21) 
(22) 
(23) 

(24) 
(25) 
(26) 



P = Po- 



(27) 
(28) 



p/f,0 

Cs = Cs,o( — )/{- ), 

tH tH,0 

where pH,o and rta,o are the background density and the 
turn-around radius at the present epoch, respectively. We 
substitute equations (27) and (28) for (11), and define Ao as 

£ ( Pa \\-Ai '^s,otH,o \2{\-B) 
. ''^PH,o' ^ rta.o ' 1-A 2(l-S),2S-3 

= 7 Ph '^ta ^H , 



7(7 - 1) 

Ko 1~A 2{l-S),2S-3 

:Ph ^ta '^H 



^(^ "H ■ (29) 

Substituting equations (21)-(26), and (29) for (17)-(20), the 
fluid equations become. 



(30) 
(31) 



(32) 

2 [(C-l)+ (^-^^) ^] + 2-(^'^^)'^^ 



7-1 



C dX' D dX 

D^'^C'^'^''^ +-in. (33) 



(7-1) 

The non-dimensional mass M is derived by another inte- 
gration(Bertschinger 1985). The equation of continuity is 
rewritten as follows: 



diT^ . 2 

— h 47rr pv — 0. 

dt 



(34) 



Substituting equations (21)-(25) for (17) and (34), the di- 
mensionless mass is 



M 



3C-2 



x^D(y -s,x). 



(35) 
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Figure 2. Three types of solutions with e=l (n=3), H=0, A=2, 5=— 7/2, and 1^0=0.01. Each hne shows the resultant profiles of the 
'stagnation' solution (dashed line), 'eigensolution' (solid line), and 'adiabatic' solution (dotted line), respectively. The shock radii are 
Xs=0.338976, 0.326595, and 0.3, respectively. Each panel represents radial profiles of the density (top left), the proper velocity (top 
right), the square of the sound speed corresponding to the temperature (bottom left), and the ratio of the Hubble time to the cooling 
time (bottom right), respectively. 



2.4 Numerical integration and boundary 
conditions 

Assuming that the gas pressure is zero, the physical variables 
of a gas shell are given by the parameter ^(Peebles 1980): 



^ _ 9 {e- sinef 



2,6',,6'-sin6' ; 

I ) ' 



v = x 



2(l-cos^)3(l + 3ex)' 
sin6l(6' - sine) 



(1 - cose)^ 



C = 0, 



2 (1- cos 61)3 



(36) 

(37) 

(38) 
(39) 
(40) 



with x=l — 3V/2X. The outer boundary condition is defined 
at the turn-around radius, 0=n. The collisional gas shell falls 
toward the center after the turn around, and the central gas 
density increases. Eventually, adiabatic compression occurs, 
and a shock propagates outward. The shock radius appears 
at a fixed fraction of the turn-around radius. The shock jump 
conditions are written as 



D2 



7 + 1 



Di, 



7-1 
7 + 1' 



27(7- 



(7 + 1)^ 



-{Vi-^XsY 



(41) 
(42) 
(43) 



where subscripts 1 and 2 represent the pre- and post-shock 
quantities, and Xs is the non-dimensional shock radius. We 
need to determine the shock position. Bertschinger(1989) 
and Abadi et al.(2000) showed three types of solutions by 
changing the shock radius. Figure 2 represents the three 
possible solutions of the 'stagnation' solution (dashed line), 
'eigensolution' (solid line), and 'adiabatic' solution (dotted 
line) for e=l(n=3), n=0, A=2, B=-7/2, and iro=0.01. 
Each shock radius is Xs=0. 338976, 0.326595, and 0.3, re- 
spectively. The ratio of the Hubble time to the cooling time 
at each radius is shown in Fig. 2, which is defined by 



^cool 



(44) 



The first type is the 'stagnation' solution (dashed 
line), which diverges at a small radius. We found that the 
similarity solution with =0.338976, which corresponds 
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to the shock position of Bertschinger's adiabatic solution 
(Bertschinger 1985), becomes the 'stagnation' solution. Wc 
analytically explain that the divergence at the small radius 
occurs due to cooling. Because the profiles of the fluid vari- 
ables arc not influenced by cooling in the outer region, we 
assume that the velocity approaches V=VoX, which is the 
same as that of the adiabatic solution (Chuzhoy & Nusser 
2000). The density and temperature can be expressed from 
equations (31) and (33) as follows: 

D = DoX% C = CoX\ (45) 

where 

* = V^[-(^ - 1) - ^^"-^C-- - |(7 - l)Vo 

^ih^H]. (47) 

As the gas relaxes to a hydrostatic equilibrium in order to 
vanish the velocity as X ^ 0, 2t — s=2 is satisfied by equa- 
tion (32), in which case Vb is 

Vo = ^:^[-KoD^-'C"'-^ + 7(7 - m]- (48) 

The first and second terms on the right-hand side represent 
the cooling and heating terms, respectively. We consider the 
behavior of the cooling term in the inner region in the case 
of Tl—Q. Assuming that the cooling term converges to zero 

at the center, /)'4-i(^2s-2 ^gggmes 

where 

/t=-|(A-l)-Fi^(2B-2) >0. (50) 

However, the parameters A and B, which satisfy k >0, are 
not in the range of < e < 1. Thus, k < is satisfied and the 
first term on the right-hand side of equation (48) diverges, 
resulting in a steeper drop of the temperature. That is to 
say, a cooling catastrophe occurs in the inner region. 

The second type is an 'adiabatic' solution (dotted line). 
It is a free-fall solution to a point mass. Thus, the den- 
sity and velocity asymptotically approach DoaX"^^^ and 
V<x.X~^^'^ , respectively. Because the cooling and heating 
times axe longer than the flow time, they are negligible for 
X^Q. For the law conservation of energy, the sound speed 
approaches C(xX~^/'^ as X — »0. 

The third type of solution is the 'cigensolution' (solid 
line), which is a marginal case of 'stagnation' solutions. 
That is to say, the 'eigensolution' satisfies the condition 
tfiom='tcooi at the origin. A divergence of the physical vari- 
ables does not occur, and the flow extends to X—Q. The 
'eigensolution' is the only solution where the mass is zero at 
the origin. We only derive the 'eigensolution' below. 



3 RESULTS 

We derive the similarity solutions of various cases of {A,B). 
The first case is {A,B)={3/2,1), which is the same as Abadi 
et al.(2000). The others are focused on the free- free and line 




Figure 3. Similarity solutions witliout cooling and heating. Each 
line indicates the results of €=1 (solid line), e=2/3 (dashed line), 
and e=l/3 (dotted line), respectively. Each panel represents radial 
profiles of the density (top panel), the proper velocity (middle 
panel), and the square of the sound speed corresponding to the 
temperature (bottom panel), respectively. 

emissions. As we can see from Fig. 1, the free- free and line 
emissions are out of the range < e < 1. Thus, we fix one 
parameter to -5=1/2, —1/2, and A=2, in which case the 
other parameter, A or B, is obtained from equation (13). 

3.1 Adiabatic solution 

We first represent the solutions without cooling and heat- 
ing for a comparison. Figure 3 shows the solutions for e=l 
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e 




u 


V 


1.0 


0.338976 


-2.27 


-0.29 


2/3 


0.289976 


-2.05 


-0.09 


1/3 


0.188952 


-1.72 


+0.14 



Table 1. Properties of similarity solutions without cooling and 
heating. Xg is the non-dimensional shock radius, u and v are 
the gradient of the density and the square of the sound speed 
corresponding to the temperature which are fitted from X=0.01 
to 0.04. 



(solid line), e—2/3 (dashed line), and e=l/3 (dotted line), 
respectively. Each shock radius is Xs=0. 338976 for e=l, 
0.289976 for e=2/3, and 0.188952 for e=l/3, respectively. 
e=l/3(n=— 1) is the closest to the slope of the CDM power 
spectrum on cluster scales (Tadros et al.l998). Outside of 
the shock radius, the gas falls freely to the center. After 
passing through the shock front, the gas velocity becomes 
nearly zero, and the gas is approximately in a hydrostatic 
equilibrium. Chuzhoy & Nusser (2000) analytically obtained 
asymptotic slopes of the similarity variables, and found 
D<xX-'-^'"+^^''^"+''^ and TocX-("-i>/("+5) as X-*0 in the 
case of e>l/6(n>— 2). For e=l, 2/3, and 1/3, the density 
becomes oc^-^ ^s for e=l, ocX'^ for e=2/3, and ocX-i-^ 
for e=l/3, and the temperature becomes ocX~° '^^ for e=l, 
ocX'^ for e—2/3, and (xX"''' for e=l/3, respectively. From the 
results of a rmmcrical calculation, we confirmed that their 
slopes do not match the analytic slopes because the velocity 
is not completely zero. Table 1 gives the shock radius Xs, 
as well as the slopes of density u and temperature v, which 
were fitted from X=0.01 to 0.04. In the case of e=l/3, the 
temperature shows a decreasing profile toward the center. 
The reason for this is the strength of the shock wave. Be- 
cause the shock velocity is Va= dr/dt oct^~^, the shock wave 
becomes strong with time when ^ >1.0 (e <2/3), resulting 
in a decreasing temperature profile toward the center. 



3.2 Case of A=3/2, B=l 

We investigate the properties of solutions for A oc p^^^T, 
which is the same cooling function as that of Abadi et 
al.(2000). Figure 4 shows the resultant profiles for 7fo=0.01 
without heating, Ti.—0. Each line shows the results of e=l 
(solid line), e=2/3 (dashed line), and e=l/3 (dotted line), re- 
spectively. Table 2 gives the shock radius Xs , and the slopes 
of the similarity variables which were fitted from X=0.01 
to 0.04. Outside of the shock radius, the gas pressure is as- 
sumed to be zero. Thus, cooling and heating are neglected 
in this region. The fluid variables have the same shapes as 
the adiabatic similarity solution in Fig. 3. For X < Xs, the 
gas has almost hydrostatic distributions. However, in the 
inner region, where the radiative cooling time is less than 
the age, radiative cooling affects the gas distribution. As a 
result, the pressure support decreases, and a cooling inflow 
is established. In Fig. 4, the non-dimensional changes of the 
internal energy and the work done on a unit mass are shown, 
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n 


Xs 


u 


V 


w 






Ko-- 


=0.01 











0.311328 


-2.03 


-0.81 


-0.02 


1.0 


0.4 


0.385248 


-1.88 


-0.64 


-0.36 




0.8 


0.519181 


— 1.15 


+0.60 


+1.12 




u 




— Z.UU 


— U.D / 


1 n no 


2/3 


0.4 


0.320795 


-1.79 


-0.41 


-0.42 




0.8 


0.397242 


-1.26 


+0.50 


+1.27 




U 


U.l ( yiD4 


1 OQ 


— U.z / 


1 n 1 c: 
+U.iO 


1/3 


0.4 


0.199747 


-1.64 


+0.03 


-0.02 




0.8 


0.224469 


-1.33 


+0.38 


+2.06 






Ko 


=0.1 











0.185853 


-2.01 


-0.97 


+0.002 


1.0 


0.4 


0.209495 


-1.99 


-0.96 


+0.03 






n 9471 0^ 


1 Qfi 
— 1 .yu 




-1-0 OR 







0.154968 


-2.00 


-0.94 


-0.001 


2/3 


0.4 


0.171938 


-1.98 


-0.93 


-0.03 




0.8 


0.197048 


-1.95 


-0.92 


-0.06 







0.099695 


-1.97 


-0.88 


+0.003 


1/3 


0.4 


0.107732 


-1.94 


-0.87 


+0.02 




0.8 


0.118345 


-1.92 


-0.85 


-0.06 



Table 2. Properties of solutions for A oc p^/^T. H is the param- 
eter of heating, w is the slope of the velocity fitted from X=0.01 
to 0.04, and the other parameters are the same as those given in 
Table 1. 



which are defined by 

dE _ rj, 
dt tj, ' 



-Pdjl/p) 
dt 



' ta 



n, 



(51) 



(52) 



where P is the gas pressure. Using equations (21)-(24), the 
non-dimensional variables are given by 



9: 



2C^ 



7(7 - 1) 



V -jX dC 
C dX^ 



(53) 



n = — [ 



■2]. 



(54) 



C V -jX dP 
7 ^ L> dX 

In Fig. 4, the gas is compressed because 11 >0. More than 
60 percent of the work is turned into internal energy, and 
the gas heats up. Since G increases with decreasing radius, 
the temperature increases with decreasing radius. 

We discuss the slopes of the similarity variables, includ- 
ing cooling. In Fig. 4, the density slopes of e=l and 2/3 
are shallower than those of Fig. 3. For e=l/3, the density is 
steeper than the result of Fig. 3. On the other hand, the tem- 
perature increases with decreasing radius in all cases. Abadi 
et al.(2000) found that the slopes of the density and temper- 
ature approach —2 and —1 as A"^0, respectively. In the case 
of {A,B), assuming that the infall velocity is V «—(,X and 
approach V=VoX^ as X— >0, the density and sound speed 
are expressed from equations (31)-(33), and (35) as follows: 
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Figure 4. Similarity variables with A=3/2, B=l, i<'o=0.01, and 1-1=0. Each line shows the solutions of e=l (solid line), e=2/3 (dashed 
line), and e=l/3 (dotted line), respectively. Each panel indicates radial profiles of the density (top left), the proper velocity (top right), 
the square of the sound speed corresponding to the temperature (middle left), the ratio of the Hubble time to the cooling time (middle 
right), the change of internal energy on the unit mass (bottom left), and the work done on the unit mass (bottom right), respectively. 



D = DoX% 
where 

w + B-2 



A-1 
B-4: 



C = CoX\ 



(55) 



(56) 



t = (57) 
From Table 2, in the case of Ko=0-l, the infall velocity is ap- 

© 2002 RAS, MNRAS 000, 1-12 



proximately constant (i.e.,«;=0). Substituting A=3/2, B=l, 
and w=0 for equation (56), the density slope becomes —2, 
and the temperature slope is —1 from equation (57). We can 
see form Table 2 that the density and temperature slopes of 
numerical calculations for Ko=0.1 are considerably close to 
—2 and —1, respectively. In the case of A"o=0.01 without 
heating, the infall velocity is approximately constant. How- 
ever, V«—^X is not fully satisfied because the energy loss 
is smaller than that of Ko=0.1, in which case the gas distri- 
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logio(X) logio(X) 

Figure 5. Similarity variables with ^4=3/2, B=l, _ftro=0.01, and ^=0.8. The middle right panel indicates the ratio of the cooling rate 
to the heating rate, and the other panels are the same as in Fig. 4. 



bution approaches the adiabatic one. Thus, the density and 
temperature slopes turn into the middle gradient between 
the case of Ka~0.1 and the adiabatic solution in Fig. 3. 

Next, we represent the similarity solutions that incor- 
porate the effects of heating and cooUng. Figure 5 shows the 
resultant profiles for ?i=0.8, and the other parameters have 
the same values as in the previous model. We show the ratio 
$ of the cooling rate to the heating rate, 

A ^ Ku 

r H7(7-l)' 



£jA-1^2S-2 



(58) 



Since the heating rate is greater than the cooling rate be- 
hind the shock, the infall velocity is weakened by heating. 



and the gas expands because 11 <0 and cools down as Q <0. 
The slope of the density becomes substantially shallow com- 
pared with the result of no heating, shown in Fig. 4, because 
the concentration of mass is suppressed due to heating. On 
the other hand, temperature increases with decreasing ra- 
dius for X>0.1. However, the change of internal energy (G) 
steeply decreases with decreasing radius, except in the cen- 
tral region. As a result, in X=0.003 ^ 0.1, the temperature 
becomes decreasing profiles toward the center. When the 
flow approaches the central region, the cooling rate becomes 
greater than the heating rate. In this region, radiative cool- 
ing produces a cooling inflow, the velocity becomes negative 



© 2002 RAS, MNRAS 000, 1-12 



Self-similar collapse with 



e 


A,B 


n 




u 


V 


w 


Ko=0.01 
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+0.34 
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2/3 


3/2,1/2 





0.277630 


-2.11 


-0.44 


+0.24 
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0.397755 


-1.26 


+0.46 


+1.64 


1/3 


4/3,1/2 





0.186403 


-1.79 


+0.08 


+1.06 






0.8 


0.229463 


-1.11 


+0.44 


+0.86 


1.0 


5/3,-1/2 





0.323424 


-2.31 


-0.54 


+0.38 






0.8 


0.496057 


-1.71 


-0.07 


-1.36 


2/3 


3/2,-1/2 





0.283470 


-2.12 


-0.23 


+0.55 






0.8 


0.398342 


-1.25 


+0.40 


+2.61 


1/3 


1.0,-1/2 





0.188644 


-1.73 


+0.14 


+1.47 






0.8 


0.232110 


-1.00 


+0.42 


+0.94 


Js:o=o.i 


1.0 


14/9,1/2 





0.230693 


-2.23 


-0.91 


+0.22 
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0.307894 


-2.14 


-0.86 


+0.10 


2/3 


3/2,1/2 





0.210950 


-2.25 


-0.82 


+0.29 






0.8 


0.273013 


-2.14 


-0.74 


+0.13 


1/3 


4/3,1/2 





0.166035 


-2.10 


-0.30 


+0.56 






0.8 


0.201210 


-1.76 


-0.02 


+0.13 


1.0 


5/3,-1/2 





0.262117 


-2.48 


-0.77 


+0.54 






0.8 


0.334214 


-2.35 


-0.71 


+0.32 


2/3 


3/2,-1/2 





0.247433 


-2.38 


-0.52 


+0.65 






0.8 


0.313539 


-2.16 


-0.36 


+0.15 


1/3 


1.0,-1/2 





0.185988 


-1.76 


+0.11 


+1.36 






0.8 


0.226802 


-1.11 


+0.42 


+0.88 



Table 3. Properties of solutions which fixed the temperature 
dependence of the radiative cooling functions to B=l/2 or —1/2. 
The other parameters are the same as those given in Table 2. 



and the gas is compressed. Because the change of the in- 
ternal energy (G) shows a steeper rise toward the center 
in Jf<0.003, the gas heats up and the temperature shows 
increasing profiles toward the origin. 

For the case of Ko~0.1, wc can sec from Tabic 2 that the 
slopes of the density and the temperature do not change sig- 
nificantly by increasing the coefficient of heating, Ti, and are 
about —2 and —1, respectively. The reason is that the flows 
are not affected by heating because the cooling is stronger 
than heating at all radii. 

3.3 Cases of 5=1/2, or -1/2 

We derive the similarity solutions fixed to B=l/2 or —1/2. 
The former is similar to the temperature dependence of 
free- free emission. Since A is determined by equation (13), 
A=U/9 for e=l, ^=3/2 for e=2/3, and ^=4/3 for e=l/3, 
respectively. The latter is similar to the temperature de- 
pendence of the line emission. Similarly, A=5/3 for e=l, 
A=3/2 for e=2/3, and A=l for e=l/3, respectively. Because 
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logio(X) 



Figure 6. Similarity variables with e=l/3(n=— 1), A=4/3, 
B=l/2, and ^0=0.01. The solid and dashed lines indicate the 
results without heating, 7i=0, and with heating, 7^=0.8, respec- 
tively. Each panel represents radial profiles of the density (top 
panel), the proper velocity (middle panel), and the square of the 
sound speed corresponding to the temperature (bottom panel), 
respectively. 



the realistic free-free and line emissions are proportional to 
the square of the density, A=2, the energy loss of the cool- 
ing functions assumed here is less than that of the realistic 
one. Figures 6 and 7 show the solutions of {A,B)— (4/3,1/2) 
and (1,-1/2) with e=l/3(n= — 1) and _R'o=0.01, respectively. 
The solid and dashed lines display the solutions without 
heating, Ti=0, and with heating, 7^=0.8, respectively. The 
shock radius Xs , and the slopes of the similarity variables 
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Figure 7. Same as in Fig. 6, except A=l and B=—l/2. 



arc given in Table 3. For e=l/3, the energy loss of cool- 
ing is small because ^=3/4 and 1. Thus, the density and 
temperature without heating are almost the same profiles 
as the adiabatic solutions in Fig. 3. For the heating solution 
with 7Y=0.8, the slope of the density becomes shallow and 
the temperature is lower than that of no heating, Ti.=0. The 
main reason is the same as the previous discussion in Section 
3.2. 

From Table 3, for e=l and 2/3 in the case of Ko=0.1, 
the density slopes are steeper than those of K(xp^^^T. If the 
infall velocity is V «—£,X, the slope of temperature ap- 
proaches —1, from equation (57). On the other hand, the 
density slope depends on the values of j4, B, and the slope 



Figure 8. Similarity variables with A=2, B=5/2, and Ko=10^'^ ■ 
The dashed line indicates the result with heating, 7^=0.8. The 
panels indicate the same variable as in Fig. 6. 



of velocity w from equation(56). For example, substituting 
(yl,B,w)=(14/9,l/2,0.22) and (5/3,-1/2,0.54), which is the 
case of £=1 with 7^0=0.1, for equation (56), the density 
slopes become —2.3 and —2.94, respectively. The resultant 
slope of the numerical calculation is shallower than this slope 
because V«—^X is not fully satisfied. 



3.4 Two-particle processes, A— 2 

We consider the similarity solutions for the case that cooling 
function has a coUisional process, A=2. We exclude the case 
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0.4 
0.8 


0.378087 
0.438560 


—2.32 
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-0.70 
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Ko- 


=0.1 






1.0 




0.4 
0.8 


0.283055 
0.304111 
0.323845 


-2.48 
-2.36 
-2.06 


-0.52 
-0.47 
-0.51 


+0.91 
+0.59 
-0.06 


1 4. Same as 


Table 2, except A=2 and B=- 


-7/2. 


Ko 


n 


Xs 


u 


V 
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10-4 


0.7 
0.8 


0.224706 
0.231401 


-1.23 
-1.10 


+0.42 
+0.45 


+0.80 
+0.82 


10-5 


0.5 
0.8 


0.214030 
0.232478 


-1.32 
-0.91 


+0.36 
+0.33 


+0.52 
+0.20 


10-" 


0.3 
0.8 


0.203330 
0.232689 


-1.49 
-0.96 


+0.28 
+0.39 


+6.24 
+1.18 



Table 5. Properties of solutions for e=l/3(n=— 1), A=2, and 
B=5/2. The parameters are the same as those given in Table 2. 



of e=2/3 because of A=3/2 from equation (13). Thus, we 
derive the similarity solutions with e=l and 1/3. For e=l, 
B=—7/2. Tabic 4 gives the results of the slopes of the sim- 
ilarity variables. The density profiles are steeper than those 
of Aocp^/^T. Substituting {A,B,w)={2 -7/2,0.91), which is 
the case of Ko=0.1 and 7^=0, for equation (56), the density 
slope becomes —4.59, and the temperature slope becomes 
— 1 from equation (57). However, because B=—7/2, the en- 
ergy loss is not large enough to produce a large infall veloc- 
ity. As a result, the gradients of the density and temperature 
become about —2.5 and —0.5 in the case of Ko=0.1. 

For €=l/3(n=-l), the cooling function is Acx p^T^^'^, 
in which case we only find the 'stagnation' solution in the 
cases of Ko=Q.Ql or 0.1. The reason is that the effect of cool- 
ing is very strong at the central region, leading to a cooling 
catastrophe. Therefore, we try to find heating solutions that 
can balance the cooling. As a result, when Ko = lQ^'^ , we find 
that the heating solutions exist with 7^=0.7-0.8, and that the 
heating is so weak that a cooling catastrophe is unavoidable 
in the case of TL <0.7. Similarly, when K{)=10~^ , an 'eigen- 
solution' exists with 7i=0.5-0.8. When _R'o=10^'', we find an 
'eigensolution' with 7^=0.3-0.8. Table 5 gives its properties. 
When Ko >1Q~'^ and Ti. >0.8, the heating is strong in the 
outer region where cooling can be neglected, and the ex- 
panding velocity becomes large and an 'eigensolution' is not 
found. Figure 8 shows the result of e=l/3, i^o=10~'', and 
7^=0.8. Because the effect of heating is strong in the outer 
region, the slope of the density becomes shallow and the 



temperature decreases with decreasing radius, which are the 
same reasons as described for the previous model discussed 
in Section 3.2. 



4 DISCUSSION AND SUMMARY 

In this paper, we present self-similar solutions including 
cooling and heating for coUisional gas. We assume that the 
cooling rate has a power-law dependences of the gas den- 
sity and temperature, Aocp^T^, and that the heating rate 
is FocpT. In order to obtain a similarity solution, A and B 
are selected by requiring that the cooling time is propor- 
tional to the dynamical time. The main results in this paper 
are the following: 

(i) In the region where the cooling rate is stronger than 
the heating rate, a cooling inflow is established, and the 
gas is compressed and heats up. Because the compression 
is stronger in the inner region than in the outer region, the 
temperature increases with decreasing radius. Furthermore, 
if the large infall velocity, V«—i^X, is produced due to 
an enormous energy loss, the temperature slope approaches 
— 1, and the density slope approaches a value that depends 
on A, B, and the velocity slope w. 

(ii) In the region where the heating rate is larger than 
the cooling rate, the infall velocity is suppressed by heating, 
the compression is weakened and the gas cools down. The 
slope of the density becomes shallow due to suppression of 
the collapse, and the temperature is lower than that without 
heating. 

(in) For e=l/3 (n=-l) and (A,B)=(2,5/2), the 'stagna- 
tion' solutions are only derived in the cases of 7^0=0.1 or 
0.01. The reason for this is that the cooling is very strong 
in the inner region, and a cooling catastrophe occurs. We 
include the heating and a small Ko in order to decrease the 
effect of cooling, and find heating solutions that can balance 
the cooling with 7!'o=10"'' and 7^=0.7-0.8. 

We consider that the self-similarity solution presented 
here can be used for modeling of structure formations. For 

example, the value of the slope parameter, 13 fu, obtained 
from observations of the surface brightness profile indicates 
the typical value for rich clusters, Pfu ~ 2/3. On the other 
hand, group systems have significantly fiatter slopes with 
fifit ~ 0.4 (Arnaud & Evrard 1999; Helsdon & Ponman 
2000). As an origin of such discrepancies, it has been pro- 
posed that the intracluster medium (ICM) extends to the 
outer regions by non-gravitational heating, leading to flat- 
ter density profiles (Metzler & Evrard 1997). Furthermore, 
in the model A=2 and e=l/3(r7,= — 1), which is close to those 
of real clusters, a cooling catastrophe occurs if there is no 
heating, and there are the heating solutions that can bal- 
ance the cooling. This may mean that real clusters oscillate 
between a heating dominated regime and one with a cool- 
ing catastrophe (Kaiser & Binney 2003) , or that the heating 
sources, such as AGN (Churazov et al. 2001; Quilis et al. 
2001; Reynolds et al. 2002; Briiggen et al. 2002; Briiggen & 
Kaiser 2002; Briiggen 2003) or thermal conduction (Fabian 
et al. 2002; Ruszkowski & Begelman 2002; Voigt et al. 2002; 
Zakamska & Narayan 2003), balance the cooling. Of course, 
it is diSicult to directly apply the solution presented here to 
a real cluster, because the time evolutions of the cooling and 
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heating differ from the Hubble time. For a realistic treatment 
of this problem, it is necessary to perform a numerical simu- 
lation while assuming regular radiative cooling and heating 
rates. However, the accuracy of the numerical simulation is 
limited by the finite resolution. Wo consider that a compari- 
son between the simulation and the similarity solution gives 
insights into the thermal evolution of clusters of galaxies. 
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